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I. INTRODUCTION 


Solar cell arrays constitute the major source of reliable long-term 

power for spacecraft orbiting the earth. High voltage solar arrays are 

desired to optimize both spacecraft mass and power efficiency. The 

1-3 

space plasma environment, though, can result in large currents being 

collected by exposed solar cells, with corresponding reductions in power 

output from the array. The obvious solution of using a covering of 

transparent insulation is at least partially offset by the expectation 

of defects, either from the manufacturing process or resulting from 

collisions with micrometeors. Early experiments showed that positive 

electrodes behind pinhole openings in the insulating sheets could 

collect electron currents far in excess of what would be expected from 

4 

electrostatic theory. Subsequent investigators verified these large 

currents, but showed a wide range of results, depending on materials, 

5—8 

configuration, and operating conditions used in the tests. 

The purpose of this research is to investigate the phenomenon of 
unexpectedly large leakage currents collected by small exposed areas of 
high voltage solar arrays operating in a plasma environment used to 

simulate space. This report covers the progress since the last annual 

9 9-10 

report. The results of the two previous annual reports and some of 

the material for this support period is reviewed in the paper of 

Appendix 

The experimental research done during the present grant period 

involved only one type of insulating material, polyimide (Kapton) . 

Positive bias tests (relative to the plasma) were done studying the 

9 

glow discharge mode of current collection. Tests were performed 



2 


studying the transition from space-charge flow current collection to 

current collection well above space charge flow values* Variation of 

current collection with time was also investigated. 

Some modeling for electron collection was done. This involved 

some numerical solutions, but is not a detailed computer simulation of 

12 

the phenomenon as others are attempting. 

Appendix B is a summary of some of the properties of secondary 

electron emission from insulators. Several of these properties have 

9-10 

been mentioned in previous reports and some are included for 
possible use in future modeling of electron collection. 
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II. APPARATUS AND PROCEDURE 
Vacuum Facility 

The experiments were conducted in a 45 cm diameter, 72 cm long 
cylindrical vacuum system. An argon hollow cathode mounted hori- 
zontally at the base of the vacuum system was used to generate the 
plasma. The hollow cathode supplied electrons which then ionized the 
argon. The only source of argon gas in the vacuum system was the flow 
through the hollow cathode. This is the same configuration used 
previously. ^ 

The operating pressure in the vacuum system was in the range of 

9-100 X 10 ^ Torr. The electron temperature ranged from 3.2 to 7.9 eV 

5 -3 

and the plasma density ranged from 1-26 x iQ cm . The range of 

plasma parameters approximates the environments expected for both an 

electrically propelled spacecraft, due to its own charge exchange 
2 

plasma and a spacecraft in a natural low earth orbit plasma environ- 
3 

ment. Measurements of the plasma characteristics were taken with a 

spherical Langmuir probe. The measurement techniques have been described 

4 

in detail elsewhere . 


Solar Array Simulation 

The insulated array of solar cells were simulated by placing a 
sheet of insulating material (polyimide) over a conductor. The sheet of 
polyimide was held to the conductor by a pressure sensitive adhesive 
(Y966 by 3M) . A small hole was punched into the insulating sheet and 
adhesive to simulate the defect. A heater and thermocouple were 
attached to the conductor. This permitted the control and measurement 
of the conductor temperature. A more detailed description is given in 
the previous report . ^ 



Procedure 


The simulated solar cell was centrally and vertically mounted in 
the vacuum system. A plasma was generated and the temperature of the 
sample adjusted to the desired setting (room temperature, 300®K, 
unless otherwise specified) . The sample was then bombarded by the 
plasma for at least one hour to stimulate desorption of gases on the 
insulator surface. A slowly increasing voltage was applied to the 
conductor, while the current collected through the hole was monitored 
along with the conductor potential. Langmuir probe characteristics 
were taken for determination of plasma properties after each test. 
Using known correlations, the plasma properties are set to approxi- 
mately the desired values before each test. Test variables included: 
hole diameter, polyimide thickness, conductor temperature, background 
pressure, and voltage magnitude and polarity. 
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III. EXPERIMENTAL RESULTS 
R. P. Stillwell 

Positive Bias - Electron Collection 
Normalization . The current /voltage data has been normalized by 
dividing the total hole current by the random current density,^ 

(3-1) 

where e is the magnitude of the electronic charge, n is the electron 

density, k is Boltzmann's constant, T^ is the electron temperature, and 

m is the electron mass. All values are in SI units, unless otherwise 
e 

specified. This normalization was found to successfully eliminate 

scatter in the data due to electron density and temperature differences 

2 

as reported previously. 

3 

Effect of Insulator Area . In the first support period, data were 
taken using different insulator areas. These data were presented in a 
logarithmic current plot, and it was concluded that there the antici- 
pated strong effect of area was not evident. The tests used three sample 

sizes: 2x2 cm, 5x5 cm, and 12 x 12 cm. These samples gave a surface 

2 2 

area range of 4 cm to 144 cm . The results of these tests are again 

shown in Fig. 3-1. From this linear plot, it is seen that some area 

effect is evident. The effect shown in Fig. 3-1, though, is considerably 

4 

smaller than that found in an earlier investigation. The normalized 
results shown lie within the extremes of the data of the earlier 
investigation, which exhibit a much steeper rise in current collection 
with increased insulator area surrounding the hole. The difference is 
possibly due to differences in experimental conditions. The earlier 
investigation used a directed ion beam composed of nitrogen for the 




laiized hole current, I /j 



areas . 
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plasma source, an electron density of 2-17 x 10 cm and a background 
pressure range of 1-10 x 10 ^ Torr. These conditions differed from the 
environment used herein (see Section II). 

The double-valued curve in Fig, 3-1 for the 5 x 5 cm sample is not 
completely understood; although this behavior is not typical, it is 
felt to be due primarily to lack of regulation in the high voltage 
power supply. 

2 

High Current Measurements . As reported previously. In collecting 
high currents, up to 5 inA, there appears to be a change in collection 
modes, indicated by the dashed lines in Fig. 3-2. The mode change that 
occurs at 1.5 mA in Fig. 3-2 is characteristic of high current measure- 
ments. During tests a bluish glow was observed to increase in intensity 
and extent with this change in mode. From these observations, a 
possible explanation for the mode change is a glow discharge initiation 
involving background neutrals. The mode change that occurs at approxi- 
mately 0.25 mA in Fig. 3-2 is less intense than the mode change at 
1.5 mA. It also is unusual in that all other mode changes observed 
were associated with luminosity changes. The lower mode change (0.25 
mA) does not have a glow associated with it, and is as yet unexplained. 

The data taken after the first mode change (0.25 mA) are still in rough 

4 

agreement with the data taken by Kennerud, indicating that any glow 
discharge effect below the glow discharge mode is probably negligible. 

Tests were performed to determine if there is any relation between 
the glow discharge observed at high currents and the background neutral 
density. The results are presented in Table 3-1. The data of Table 
3-1 indicates that the hole power (collected current times conductor 
potential) at the mode change (occurs at 1.5 mA in Fig. 3-2) increases 
with decreasing neutral density, shown in Fig. 3-3. Below about 



Hole current , mA 
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5.0 

4.5 

4.0 

3.5 

3.0 

2.5 

2.0 



Hole diameter = I mm 
Conductor temperature =40®C 
Polyimide thickness = 0.127mm 
Electron density=6.8xl0® cm* 
Electron temperature = 4.5eV 



Conductor potential, V 


Fig. 3-2. Electron current collection characteristics extended to 
higher currents. 
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Table 3-1. Results of Neutral Density Tests Involving Mode Change Point. 


Neutral 

Density 

(cm”3) 

Conductor Potential 
at Mode Change 
(volts) 

Collected Current 
at Mode Change 
(yA) 

Electron 

Density 

(cm“3) 

12 

2.9 X 10-^^ 

850 

3000 

5.6 X 10 

12 

3.2 X 10 

770 

4200 

- 

12 

4.0 X 10 

790 

290 

- 

9.7 X 10^^ 

510 

270 

1.8 X lo' 

13 

1.6 X 10 

565 

190 

1.2 X 10 

13 

1.9 X lO-^-^ 

485 

110 

2.2 X 10 

13 

2.4 X lO-^-^ 

445 

145 

1.5 X 10 

13 

3.2 X 10 

480 

125 

2.6 X 10 



Hole power at mode change, watts 
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Fig. 3-3. Hole power at mode change versus neutral density. 
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12 3 

3 X 10 neutrals/cm , the required power rises to high values. The 

glow discharge process would then not.be expected to occur in space 

where much lower neutral densities are found. Future tests should then 

be performed in an environment with a significantly lower neutral 

12 3 

density than is presently used, preferably below 3 x 10 /cm . 

Effect of Hole Size at High Current Measurements . Studies were 

done to determine if current collection at high currents was affected by 

hole size. Five hole diameters were tested; 1.0 mm, 2.0 mm, 3.0 mm, 4.1 

mm, and 5.0 mm. The results of these tests are presented in Table 3-2. 

The current/voltage curves did not correlate well when using the random 

current density factor, j^, nor did they correlate well when using the 

factor of random current times hole area (I = i x hole area) . How- 

^ 0*^0 

ever, if the conductor potential at a current of 1.5 mA is evaluated, it 
is seen that the potential at this current decreases with increasing 
hole size. This result suggests that hole size is an important param- 
eter for collection at high currents in the absence of a glow discharge 
collection mode. The conductor potential at 5 mA is shown in Table 3-2. 
(The 5 mA data corresponded to a glow discharge in all cases.) No 
correlation was found among hole size, electron density and current 
collection in the glow discharge mode. 

These tests indicate that further studies of hole size versus high 
current collection, in the absence of a glow discharge collection mode, 
are necessary. 

Multiple Holes . During tests, when the laboratory was darkened, 
small localized bluish glows were observed around the holes prior to the 
mode change (glow discharge mode). To study these glows further, 
experiments were done in which two 1 mm diameter holes were punched in 
the polyimide. Two tests were performed in which the holes were 
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Table 3-2. Results of Various Hole Diameters at High Currents. 


Hole 

Diameter 

(mm) 

Electron 

Density 

Electron 

Temperature 

(eV) 

Neutral 

Density 

(cm“3) 

Conductor 
Potential 
at 1.5 mA 
(volts) 

Conductor 
Potential 
at 5 mA 
(volts) 

1.0 

1.3 X 10^ 

6.1 

12 

8.0 X lO'^'^ 

860 

440 

2.0 

5.7 X 10^ 

5.7 

12 

9.7 X 10-^^ 

720 

1050 

3.0 

1.5 X 10^ 

5.5 

12 

9.7 X 10-^^ 

500 

2200 

4.1 

1.1 X 10^ 

4.6 

12 

9.7 X 10-^^ 

335 

1900 

5.0 

1.1 X 10^ 

3.2 

12 

9.7 X 10-^^ 

320 

730 
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separated by 2 cm and 5 cm (center to center). The results show no 
interaction between holes when the glows are first observed, as indi- 
cated by measurements of total power, neutral density or electron 
density. 

The glows became visible in both holes at the same time. Other- 
wise, they appear to be independent of each other. The glows differed 
in both intensity and extent in each hole. With increasing current, the 
glows gradually grew until they merged with each other. This occurred 
at 0.33 watts (510 volts) for the 2 cm hole separation and at 1.10 watts 
(1100 volts) for the 5 cm hole separation. Once merged, the glow 
continued to expand outward. 

During the tests, arcs were also visible in the holes. The arcs in 
each hole also appeared to be independent. These arcs were not always 
visible, disappearing and returning in an apparently random fashion. 

The arcs appeared to be red at first, then becoming white. The arcs in 
the holes continued even after the glows merged. 

Negative Bias - Ion Collection 
2 

It was previously reported that for negative conductor potentials 
under a 1000-volt magnitude, collected currents agreed with values 
predicted from electrostatic theory. Tests were conducted at larger 
voltages. These tests were done with five different hole diameters; 1.0 
mm, 2.0 mm, 3.0 mm, 4.1 mm and 5.0 mm. The tests all yielded the same 
qualitative results. 

At negative potentials greater than about 1000 V, the current began 
to oscillate. These oscillations appeared to become larger with both 
time and voltage. During these oscillations, localized arcing was 
observed around the hole. As the current increased, this arcing 
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appeared more continuous and the current collection became more steady. 
The currents measured at this point were usually between 100 and 200 pA. 
Tests were stopped at about the 100-200 pA level because arcing was 
observed between the high voltage connector leading to the sample and 
the sample holder. This was typically at negative potentials greater 
than 5000 V. 

From these results there appear to be two regimes of current 
•collection. In the lower voltage regime collected currents are of the 
magnitude expected from electrostatic collection. In the higher voltage 
regime collected currents are much higher than expected from electro- 
static collection, although still considerably lower in magnitude than 
those collected with positively biased conductors. This higher collec- 
tion regime appears to be associated with arcing occurring in the hole. 
It has also been observed that a dark broimish layer is deposited on the 
conductor exposed to the plasma by the hole. This suggests that the 
polyimide is vaporized or sputtered. The transition between these two 
regimes is very unsteady and is characterized by large current oscil- 
lations and to a smaller extent by voltage oscillations. That is, when 
the arcing appears continuous, the current rises and voltage decreases. 

Variation of Current Collection with Time . Tests measuring current 
collection as a function of time, at constant negative potential, were 
performed. These tests were conducted to determine if the dark brownish 
layer deposited on the conductor affected the ion collection. 

The criterion used to determine the voltage for these tests was the 
voltage value when a hole began collecting more than 100 pA. This value 
was chosen because 100 pA is in the high current collection regime. 

Table 3-3 shows the data from which the voltage for these tests was 
determined. The tabulated values of voltage do not show a consistent 
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Table 3-3. Conductor Potentials where Ion Current Collection First 
Exceeded 100 )jA, for Various Size Holes. 


Hole Diameter 
(mm) 

Conductor Potential for Collection ^100 pA 
(volts) 

1.0 

-6100 


-3100 


-5100 

2.0 

-3800 

3.0 

-5500 


-3600 

4.1 

-4100 

5.0 

-4400 
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trend with hole diameter. The average of all values in Table 3-3 is 
about -4500 volts. This value (-4500 v) was therefore used in these 
tests. 

Three hole diameters were tested; 0.35, 3.0 and 5.0 mm. These 
three hole sizes cover the range of hole sizes used to date. 

Figure 3-4 shows the results of these tests. It was observed 
during these tests that the decrease in collected current roughly 
corresponds to the development of the brownish film. 

During these tests there was also considerable voltage/current 
oscillation (see Fig. 3-4). This indicates the need for a regulated 
power supply in future tests. 



Hole current, /tA 


17 


Pinhole 
diameter 
o 0.35mm 
□ 2.95mm 
A 4.99mm 


Voltage held constant 
at -4500 volts 



Election Election 
temperature density 
5.5eV I.OX|0®cm" 
3.9eV UX|0«cm* 
6.7eV 9.5X10® cm' 


3 

3 

3 


Fig. 3-4. Current collection as a function of time for a bias 
of -4500 volts. 



18 


IV. MODELING 

. , R. P. Stillwell 

' i I 

For positive bias (electron collection) , the approach taken has 
been that part of the insulator surrounding the exposed conductor (the 
hole) takes part in the current collection. The mechanism for the 
electrons moving across the insulator material to the hole, for collec- 
tion by the conductor, is taken to be secondary electron emission. The 
driving force for this process is a potential gradient from the hole (at 
applied conductor potential) to some radius on the insulator (near 
floating potential) which does not contribute to the electron collection. 

The process envisioned is an electron from the plasma striking the 
surface of the insulator, within the collection radius, and a secondary 
electron being emitted. This secondary electron, moving in a direction 
toward the hole, gains energy from the potential gradient and strikes the 
surface with sufficient energy to cause another secondary electron to be 
emitted. In this manner, an electron striking the insulator, within the 
collection area, appears to travel across the surface to the hole. The 
secondary electron emission coefficient would have to be unity, on the 
average, for the traveling electrons to maintain current continuity and 
conservation of charge. 

In this approach, the sheath would essentially cover only the area 
of the insulator contributing to the electron collection. This model 
is visualized in Fig. 4-1. In the figure, r^ represents the radius of 
the hole while r^ represents the total collection radius. 

Collection Area 

The first calculations done were to determine if the model was 

feasible on the basis of the necessary size of the collection radius, r . 

c 
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The current is given by the expression 



^.2 .2 eV 

^ 4 ^ 4-* kT 


9Q8 

A^ + Ai(2f) 


e 


_ 


(4-1) 


2 2 

The terra i {l+b/4+(l-b /4) eV /kT } is the collection term 
o a e 

found by Parker and Whipple for collection by a planar probe surrounded 

by a grounded conductor.^ The b in this terra is a constant, found to 

2 

be about 1.8 for the test environment used herein. is the applied 

conductor potential, k is Boltzmann^ s constant, T^ is the electron 
temperature, e is the electronic charge, and is the random current 


density, given by" 



(4-2) 


where n is the electron density and m^ the electronic mass. The term 
2 2 

{1 + b /4 + (1 - b /4) eV^/kT^} represents the geometry factor asso- 
ciated with the expanding sheath. The area is given by the term 
{A^ + (298/T)^^^^^}. is the area of exposed conductor, A^ the 

area of insulating material contributing to the electron collection, and 
T is the sample temperature (in °K) . The area A^ is multiplied by the 
factor (298/T)^^^^^ to take into account the temperature dependence cf 
the collected current (also see Appendix B) . The temperature term is 
not multiplied by the total collection area because only the collection 
by the insulator is affected by the sample temperature. The expression 
for S(V^) was the best linear least squares fit to the experimental 

4 

data (see Fig. 4-2). It was found to be 

5(v) = 2.1 + 4.3 X 10"^ V . 

a 


(4-3) 



Power of te 



Conductor potential, V 


Fig, 4-2. Power of temperature variation for electron current as a 
function of conductor potential. 
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The expression for the current (Eq, 4-1) was then matched to 14 
sets of data and the radius of insulator that contributes to collection 
was determined. The insulator radius is defined as 

r . = r - r, . (4-4) 

1 c n 

r^ was assumed to be a function of only conductor potential for 
simplicity. The calculated values of the radius, r^, were then fitted 
to a simple parabolic variation with potential. The best least-squares 
fit was obtained with (see Fig. 4-3) 

r. = 3.46 X 10 "® V ^ . (4-5) 

i a 

In the model described in the previous section, there is a 

potential gradient along the surface of the collection area. This 

2 2 

would imply that the potential in the term {l+b/4+(l-b /4) eV/kT^} 
should be replaced by a potential averaged over the collection area. 

This could reduce this term, by as much as a factor of 3, which would 
mean that r, would have to be Increased as much as by a factor of 1.7. 
Even with these considerations, these calculations still demonstrate 
that the model is feasible since the collection area, in the worst 
case, is only 50% of the sample size. 

This empirical fit was then used to calculate expected 
current/voltage curves. Typical calculated curves are shown in Figs. 

4-4 and 4-5, together with the corresponding experimental data. The 
empirical fit is shown to give reasonable agreement to experimental 
results in these comparisons, and can be used to give an order of 




Normalized hole current, I/j 
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Conductor potential, V 


Fig. 4-4. Comparison of experimental and calculated electron current 
collection characteristics for a hole diameter of 1 mm, 
electron density of 7.4 x 10^ cm”^, electron temperature 
of 7.7 eV and conductor temperature of 320°K. 



Normalized hole current, I/j 
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Conductor potential, V 


Fig. 4-5. Comparison of experimental and calculated electron current 
collection characteristics for a hole diameter of 0.52 mm, 
electron density of 2.1 x 10^ cm”^, electron temperature 
of 12.8 eV and conductor temperature of 313°K. 
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magnitude calculation or better of expected results. Because of the 
empirical fit used, a change in insulator material would require addi- 
tional experimental data. 

Electron Trajectories Across the Insulator Surface 
Further calculations were done to determine if the electron motion 
across the insulator surface (described earlier) were possible. To 
accomplish this, equations of motion for the electron moving across the 
insulator surface were needed. In the derivations of these equations, 
cylindrical coordinates (r, 0, z) are used. The radius along the surface 
of the insulator is r and z is the height above the surface of the sample. 
Surface Charge Density . The first step is to derive an expression 

for the surface charge density. It is assumed that a constant current 

2 2 
density, j in amps/m , is striking the collection area, irr^ . It is 

further assumed that the collection radius is a function of the applied 

conductor potential only. In the model, all electrons striking within 

the radius r^ will be collected, so that a line current density, j^^ in 

amps/m, can be calculated for a specific radius, r. 


j 


L 
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The surface charge density is then given by the expression 


a(r) = j^/v^ 


(4-7) 


where v 

s 


is a surface velocity dependent only on the applied conductor 


potential. Thus 
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S 

By introducing a dimensionless quality, R = normalized radius, 

the surface charge density takes the form 


a(R) 



[l-R^] 


(4-9) 


in which the term jr^/2v^ is a geometric factor dependent on the 
applied conductor potential and current density, j. In this form, a 
normalized surface charge density, 2v^/jr^ cj(R) can be plotted against 
normalized radius, R, to give a profile that is independent of current 
density or applied conductor potential (see Fig. 4-6). 

Surface Potential . With an expression for the surface charge 
density, a surface potential can be calculated across the radius r^. 
The sample potential is given by the expression^ 


2tt 


o o 


g(r' )r' dr ' 

[rVr*^-2rr’cos0’ 


(4-10) 


where is the permittivity of free space. 

In terms of the normalized radius, the expression becomes 


V(R) 



Stte V 

o s 
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(l-R*^)dR* 

’ ^-2RR' COS0 



(4-11) 


where R* = rVr and R, = r, /r . From this a normalized potential, 
c n h c 

2 

87T£ V /ir V(R), can be calculated. This was desirable since the 
integral (Eq. 4-11) could not be solved analytically, the numerical 
solution (see Fig. 4-7) would be independent of conductor potential 






Normalized surface 
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or current density. The numerical solution is shown as a function of 
normalized radius. 

To get estimates of the values of j and r^, as a function of applied 
conductor potential, the equations developed in the previous section 
(collection area) for the current (Eq. 4-1) and the empirical expression 
for r^ (Eq. 4-5) were used. The approximation ^ used to 

estimate r^ (the approximation is reasonable inasmuch as r^ << r^) . 

The current density was found using 


I = 



(4-11) 


The calculated value for the current was used rather than actual data 
because the expression for the current (Eq. 4-1) is representative of 
a large set of data. 

Once the values of j and r^ were known, v^ was calculated by 
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Jirr^ 

V(R^)-V(1) 

8ire 
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(4-12) 


^ V(R^)-V(1) 

V = — :: . 

s Q 2 V 

OTT e a 

o 

These values are tabulated for the case in which the hole radius is 

5 3 

0.5 mm, electron temperature is 7.7 eV, electron density is 7.4 x 10 /cm , 
and conductor temperature is 320®K (see Table 4-1) . 

Electric Field Along the Surface . The normalized surface potential 


was approximated by the curve 
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Table 4-1. Collection Radius, Collected Current, Current Density, 
Surface Velocity and Normalized Hole Radius for Various 
Values of Conductor Potential. 


Conductor 

Collection 

Collected 

Current 

Surface 

Normalized 

Potential 

Radius 

Current 

Density 

Velocity 

Hole Radius 

(volts) 

(m) 

(amps) 

(amps/m2) 

(m/sec) 

(V^c^ 

300 

3.11x10"^ 

1.21x10"^ 

0.398 

7.00x10^ 

0.161 

350 

4.24x10"^ 

2.52x10"^ 

0.446 

1.43x10^ 

0.118 

380 

5. 00x10"^ 

3.74x10"^ 

0.476 

2.11x10^ 

0.100 

400 

5.54x10"^ 

4.74x10"^ 

0.492 

2.64x10^ 

0.090 

425 

6.25x10"^ 

6.33x10'^ 

0.516 

3.45x10^ 

0.080 

454 

7.14x10"^ 

8.65x10“^ 

0.540 

4.61x10^ 

0.070 

491 

8.33x10'^ 

1.26x10"''^ 

0.578 

6.51x10^ 

0.060 

538 

1. 00x10"^ 

1.91x10"^ 

0.608 

9.48x10^ 

0.050 

601 

1.25x10"^ 

3.23x10"^ 

0.658 

1.52x10^ 

0.040 

650 

1.46x10"^ 

4.63x10"^ 

0.691 

2.09x10^ 

0.034 

695 

1.67x10"^ 

6.36x10"'^ 

0.726 

2.74x10^ 

0.030 

750 

1. 95x10"^ 

9.11x10"^ 

0.763 

3.74x10^ 

0.026 

800 

2.21x10"^ 

1.22x10"^ 

0.795 

4.81x10^ 

0.023 

850 

2.50x10“^ 

1.61x10"^ 

0.819 

O.OOxlo'^ 

0.020 
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Sue V 
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V(R) 


6,263 


0.038 

r3/2 


4.665 . 


(4-13) 


The electric field along the surface could then be found by 
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E (r) 
r 


3V(r) 

3r • 


(4-14) 


This yielded, in terms of normalized radius 


E (R) 
r 



8tT£ V 
o s 


3.131 


0. 0564 

• 


(4-15) 


A plot of the normalized electric field along the surface, 

Sue V /ir E (R) versus normalized radius is shown in Fig. 4-8. 

Sample Calculations . The surface charge density, surface potential 
and electric field along the surface were calculated as a function of 
normalized radius for an applied conductor potential of 601 volts. This 
was done to illustrate the magnitude of these quantities. The equations 
of motion that will be solved will be for this condition. The quantities 
are plotted in Fig. 4-9. 

Electric Field Above the Surface . To write the equations of 
motion, it is necessary to estimate the electric field above the surface 
of the sample. Two approaches are possible. One, that the electrons 
traversing the surface are influenced only by the local surface charge 
density. In this case, the electric field, E^, may be approximated by^ 


E 

z 


a/e 


(5-16) 




Electric field Surface Surface 

along surface, potential, charge density 

V/m V Coul/m* 
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The (a) surface charge density, (b) surface potential, and 
(c) electric field along the surface plotted against 
normalized radius for an applied conductor potential of 
601 volts. 


Fig. 4-9, 




35 


The alternate approach is that the electrons moving across the surface 
are affected by the sheath between the surface and the plasma. In 
this case, the potential is approximated by Child's law^ 



where z is zero on the surface of the insulator and z = z at (|) = 0, 

max 

and V(R) is given by Eq. (4-13). The electric field for Child’s law is 
then 



/ 



1 1/3 


(4-19) 


Both these fields are plotted as a function of normalized radius for a 
conductor potential of 601 volts on Fig. 4-10 (Child’s law electric 
field for z = 0) . The electric fields differ by two orders of magnitude 
and should yield very different trajectories for the electrons crossing 
the insulator surface. It should then be possible to choose the correct 
approach to the z component of the electric field after examining the 
trajectories that each approach yields. 

Equations of Motion . The equation of motion for r is given by 


,2 -eE 
dr r 


-ejr 


dt 


m 


Sttc m V 
o e s 


3.131 


0.0564 




(4-20) 


Introducing the dimensionless time variable 


T 






Sue m V 
o e o 


(4-21) 



Z component of electric field, V/m 
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and the normalized radius, R, the equation of motion for r becomes 


d^R 

dx^ 


3.131 0.0564 


R 


3/2 


R' 


5/2 


For the equation of motion for z 


^2 -eE 
d z z 


dt 
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The surface charge density electric field yields 

2 


dt^ 
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Introducing the dimensionless variable 


the equation of motion, in dimensionless form, becomes 


d^Z 

dr^ 


= -47T 
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The Child’s law electric field gives 
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There are two sets of two coupled second order differential equations 
to be solved. They were solved numerically using a Runge-Kutta fifth 
order method for a system of first order coupled differential equations. 
This was done by introducing the variables 


^ dx 


(4-28) 
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5 = ^ (^- 29 ) 

which correspond to dimensionless velocities in the r and z directions, 
respectively. The conditions solved for were for a conductor potential 
of 601 volts, an initial position on the surface of the insulator at 
the collection radius, r^, and an initial velocity in the z direction 
equivalent to 2 eV. These conditions are for an electron striking the 
insulator at r = r^ and emitting a secondary electron normal to the 
surface. An energy of 2 eV was picked to correspond to the probable 
energy of a secondary electron (see Appendix B) . 

The trajectories given by the differential equations are shown in 
Fig. 4-11. For calculated from the surface charge density, the 
trajectory took 44 hops across the insulator, with all the hops, after 
the initial hop, gaining, from the potential gradient across the surface, 
approximately 10 eV or less. These conditions would not permit the 
proposed secondary electron emission mechanism to take place. 

The trajectory yielded by the Child^s law approach to is more 
consistent with the proposed model. The energy gained by the first hop 
was approximately 60 eV, and for any Initial position for R < 1, no more 
than two hops were ever necessary to reach the hole (from further cal- 
culations, not shown). By examining the electric fields, the condition 
necessary for trajectories feasible with the proposed model is 

E » E 
r z 

Calculating the Child* s law sheath thickness, the sheath would appear to 
be a large spike on the surface of the insulator. This sheath is many 



Normalized height above surface, 



U) 

VO 



Fig. 4-11, Electron trajectories for an electron moving across the insulator for E 2 calculated 
from (a) the surface charge density and (b) Child’s law, for a conductor potential 
of 601 volts and Initial velocity equivalent to 2 eV normal to the surface. 
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times the collection radius (see Fig. 4-12, Z = z /r ). This is 
in contradiction to the proposed model. The sheath is unphysical and 
fails to provide a uniform current density to the collection area. 

These calculations show that the proposed approach has been 
incorrect in assuming that the collection area is large. The two pos- 
sible solutions within the framework of the model contradict it. The 
local charge density solution has the correct sheath (planar) , but the 
trajectories are inconsistent with the secondary electron mechanism of 
surface conduction. The Child’s law solution has the correct trajectory 
but the sheath is unphysical. This says that the correct trajectories 
cannot be obtained with a planar type sheath. The alternative to this 
is a large sheath surrounding a small collection area (see Fig. 4-13). 

A spherical sheath much larger than the collection area would allow the 
condition for the correct trajectories, » E^, to take on the much 
reduced collection area. 

The conclusion of this analysis, then, is that the sheath should be 
approximately hemispherical. The radius for secondary electron emission 
(collection radius) is much smaller than the sheath radius. Between the 
collection radius and the sheath radius, the electron trajectories 
should be nearly radial to avoid collisions with the insulator surface 
and a resultant buildup of charge on that surface. Between the collec- 
tion radius and hole radius, surface conduction takes place with the 
secondary electron mechanism. This approximate model will be studied 


and tested in future work. 
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V. CONCLUDING REMARKS 

A glow discharge electron collection mode appears at high currents 

(>1.5 inA) . This collection mode appears to be related to the background 

neutral density and as such does not simulate a process that would occur 

in a space environment. Data taken before this mode change occurs are 

in rough agreement with an earlier investigation,^ indicating that most 

of the tests taken to date are still valid for space simulation. 

Future tests should be done in an environment that has a lower neutral 

12 

density than is presently used. Data indicate a density below 3 x 10 

-3 

cm should be adequate to avoid the glow discharge mode with argon. 

Negative bias (ion collection) tests indicate that electrostatic 
collection occurs for low voltages (<-1000 volts) but another mechanism 
occurs at high voltages. This other mechanism results in a brownish 
film being deposited on the exposed conductor area, decreasing ion 
current collection. The deposited film suggests that sputtering or 
vaporization of the polyimide is taking place. The tests in the high 
voltage regime exhibited oscillation of current and voltage, indicating 
that a regulated power supply should be used in future experiments. 

The modeling of electron collection was attempted. The model 
assumed a large collection region over the insulator with a sheath 
covering approximately the same area. Calculations show the model to be 
inconsistent with the assumptions, indicating that a model with a small 
collection region and a large hemispherical sheath should be considered. 
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EXPERIMENTAL SIMULATION OF SPACE PLASMA INTERACTIONS WITH HIGH VOLTAGE SOLAR ARRAYS 

R. P. Stillwell, H. R. Kaufman, and R. S. Robinson 
Colorado State University 
Fort Collins, Colorado 80523 


Abstract 


Operating high voltage solar arrays in the 
space environment can result in anomalously 
large currents being collected through small 
insulation defects. Tests of simulated defects 
have been conducted in a 45 cm vacuum chamber 
with plasma densities of 10^ to 10^ cm“^. 

Plasmas were generated using an argon hollow 
cathode. The solar array elements were simulated 
by placing a thin sheet of polyimlde (Kapton) 
insulation with a small hole in it over a con- 
ductor. Parameters tested were: hole size, 

adhesive, surface roughening, sample temperature, 
insulator thickness, insulator area. These 
results are discussed along with some preliminary 
empirical correlations. 

Introduction 


Solar cell arrays constitute the major 
source of reliable long-term power for space- 
craft orbiting the earth. High voltage arrays 
are desired to optimize both spacecraft mass and 
power efficiency. The space plasma environment, 
though, can result in large currents being col- 
lected by exposed solar cells, with corresponding 
reductions in power output from the array. The 
obvious solution of using a covering of trans- 
parent insulation is at least partially offset 
by the expectation of defects, either from the 
manufacturing process or resulting from collis- 
ions with micrometeors. Early experiments 
showed that positive electrodes behind pinhole 
openings in the insulating sheets could collect 
electron currents far in excess of what would be 
expected from electrostatic theory. At high 
enough plasma densities and electron currents, 
the current collection process degrades into a 
visible arc. At low plasma densities, the large 
electron currents are collected in the absence 
of a visible arc, and appear to be associated 
with an essentially nondestructive surface 
phenomenon. The maximum current at high voltage 
appears, from previous investigations, to depend 
primarily on the plasma density and the area of 
the insulating surface surrounding the pinhole. 

The major emphasis of this research is the 
understanding of the apparent surface conduction 
by which the area of insulation surrounding an 
exposed conductor enhances the current collection. 
Experimental tests of high voltage current col- 
lection simulating both insulator defects and 
plasma environments are described. 


except for the few tests conducted to evaluate 
collection characteristics in the presence of 
lower quality adhesives. A small hole was 
punched or drilled into the insulating sheet to 
simulate the defect. The sample was mounted in 
the center of the vacuum system. 


The operating pressure in the vacuum system 
was in the range of 9-30 x 10”^ Torr. The 
plasma density ranged from 1-25 ?< 10^ cm“^. 

This plasma range approximates the environments 
expected for both an electrically propelled 
spacecraft, due to its own charge exchange 
plasma^ and a spacecraft in a natural low earth 
orbit plasma environment . ^ 


The procedure for making a test was to 
allow the sample to be bombarded by the plasma 
for at least one hour to stimulate desorption of 
gases on the insulator surface. The current 
collected through the hole was monitored as the 
potential on the conductor was increased. 


Results and Discussion 


Positive Bias 


High Current Measurements. In collecting 
high currents, up to 5 mA, there appears to be a 
change in collection modes. Indicated by the 
dashed line In Fig. 1. During tests a bluish 
glow was observed to Increase in intensity and 
extent with the change in mode (toward increasing 
current). From these observations, a possible 
explanation for the mode change is a glow dis- 
charge Initiation involving background neutrals. 
Because this glow discharge process is probably 
associated with the high background pressure 
(higher than in a space environment) , most of 
the data presented herein will be at a current 
below any major mode change, unless otherwise 
indicated. 

Pinhole Size. Four hole diameters were 
investigated: 0.35, 0.52, 1.0, and 2.0 mm. 

This range of diameters spanned an area range of 
over 30; 1. To isolate the effect of the area of 
exposed conductor from plasma density effects, 
the current was normalized. The random electron 
current density in a plasma is given by 8 



Apparatus and Procedure 

All tests were conducted in a 45 cm diameter, 
45 cm long cylindrical vacuum system. An argon 
hollow cathode mounted horizontally at the base 
of the vacuum system was used as the source of 
the plasma. The hollow cathode supplied electrons 
which then ionized the argon. A protected array 
of solar cells was simulated by placing a sheet 
of insulating material over a conductor. The 
Insulating sheet of polyiraide (Kapton) was held 
to the conductor by a space-qualified adhesive. 


where e is the electronic charge, n is the 
electron density, k is Boltzmann's constant, T^ 
is the electron temperature, and mg is the 
electron mass. 

It would be expected that the current 
collected through a hole in an insulating sheet 
would vary with plasma density, proportional to 
the random electron current density given above. 
Another possible normalizing relationship con- 
sidered is this current density times the area 
of the hole in the Insulating sheet, Iq = jo^hole* 


Reieiised to AIAA to publish in ill forms. 
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Fig. 1. Electron current collection 
characteristics extended to 
higher currents. 


The data were normalized using both the random 
current density and the random current, see Fig. 
2(a) and 2(b). The better correlation obtained 
with Jq Indicates that electron current collec- 
tion at these densities Is insensitive to hole 
area. 

As seen In Fig. 2(a), anomalously low cur- 
rents were initially collected with the large 2 ram 
hole; however, at higher collection currents, 
some mode change occurred that brought the data 
for this hole diameter more In agreement with 
the other data. This mode change was not assoc- 
iated with any visible glow discharge and was 
well below the other mode changes associated 
with high collection currents. It was also a 
repeatable phenomenon. At this point it can 
neither be explained nor Ignored. 

These results are not in agreement with 
those reported earlier by Kennerud,^ who found 
an effect due to hole area; however, the magni- 
tudes of the normalized current collection data 
are in rough agreement. The difference is 
possibly due to differences in experimental 
conditions. This investigation used an argon 
hollow cathode for a plasma source, had an 
electron density range of 1-25 x lO^ cm“^ and a 
background pressure range of 9-30 x 10”^ Torr. 

The previous Investigation used a directed Ion 
beam composed of nitrogen for the plasma source, 
an electron density of 2-17 x iq 3 cm“^ and a 
background pressure range of 1-10 x 10"5 Torr. 



(a) 



Fig, 2. Electron current collection charac- 
teristics for a range of hole sizes, 
normalized by (a) (b) j^,. 
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Adhesive. Althougti most of the tests were 
made with a space qualified adhesive (Y966, by 
3M) , some tests were also conducted using a low 
priced commercial brand (Scotch Double Stick, 
also by 3M) . Tests were conducted at two tem- 
peratures, 27 and 95°C. The low priced com- 
mercial adhesive was found to have modified 
sticking properties (i.e., adhesive was soft and 
pliable) after testing at higher temperatures, 
while no such change was evident for the space 
qualified adhesive. Such a change in properties 
for an organic material in a vacuum environment 
is almost always associated with outgassing. 

This outgassing should he increased at the high 
temperature. Figure 3 shows the results of 
these tests. It caii !..? seen from Fig. 3 that 
there appears to be no significant effect of the 
adhesive tests, hence probably no effect due to 
outgassing from the adhesive. Other tests with- 
out any adhesive tended to support this conclu- 
sion. Except for more data scatter, possibly 
due to varying the distance between adhesive and 
sample, the data with no adhesive agreed with 
Fig. 3. 



Fig. 3. Effect of adhesive type and sample 
temperature on electron collection. 

Fresh and Old Samples. Subsequent tests 
made on the same sample exliibited a lack of 
reproducibility compared to the initial test on 
that sample. There was a general degradation in 
current collected with repeated tests, as shown 
In Fig. 4. One possible explanation for this 
effect was that adsorbed gas on the Insulator 
surface was involved in the collection process 
and was being removed during the test. To 
evaluate this possibility, the sample was taken 
out of the vacuum chamber after several tests 
of the sample to allow adsorbed gases to collect 
on the insulator surface. When the sample was 
tested after an overnight exposure to the atmos- 
phere, results similar to those obtained with 


samples not exposed to atmosphere were found. 

The effect therefore did not appear to be due to 
adsorbed gases. 



Fig. 4. Comparison between the first, 

second and third tests of a sample. 

On visual inspection of samples after 
tests, it was observed that there was some dis- 
coloration of the insulating material about the 
hole, and the interior of the hole Itself 
appeared smoother and more polished than It did 
Immediately after drilling or punching the hole. 
It was found that previously tested holes could 
be made to increase their current collection by 
merely mechanically scraping the interior of the 
hole. This scraping was accomplished with a 
drill bit and resulted in no significant change 
in hole diameter. The effect of scraping can be 
seen in Fig, 5. A possible explanation of this 
effect Is that the interior of the hole and the 
area immediately surrounding the hole is altered 
by the energetic bombardment of electrons, and 
that scraping the Interior of the hole restores 
the hole to its approximate former state. The 
increased current collection points out that the 
condition of the hole is an Important part of 
the collection mechanism. 

Surface Roughening. Because the condition 
of the surface appears to be Important to the 
collect Ion process, the flat Insulator surfaces 
were scribed with a sharp metal stylus. This 
was done to give a better understanding of the 
effect of surface condition. Three scribed 
patterns were used: orthogonal lines (two sets 

of parallel lines at right angles to each other); 
radial lines originating from the vicinity of 
the hole; and concentric circles with the hole 
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Fig. 5. Comparison of current character- 
istics of a single hole in three 
conditions, fresh, old, and reamed. 


at the common center. As seen from Fig. 6, 
there is a significant reduction in the current 
collected by the scribed surfaces. A possible 
explanation of this effect appears to be that 
secondary electron emission is reduced by rough- 
ening a surface.^ Secondary electron emission 
has been proposed as a collection mechanism by 
previous investigations.^*^ 


. Variation of Temperature. To test for the 
effect of temperature, the samples were heated 
to five different temperatures: 27®, 47®, 70®, 

95®, and 120®C. There was a general reduction 
in collected current with Increasing temperature, 
see Fig. 7. By crossplotting the faired curves 
of Fig. 7, a temperature variation was found as 
a function of conductor potential, this is shown 
in Fig. 8. This effect seems to Indicate a 
multi-step secondary electron emission process. 
For Insulators, secondary electron emission has 
a temperature variation of T“^'^ under electron 
bombardment.^® With the probability of a single 
secondary electron emission event proportional 
to the larger negative powers in Fig. 7 

can be thought of as resulting from a number of 
secondary electron events in series. The 
increased negative powers at higher collection 
voltages further implies a larger number of 
secondary electron events in this series. 


Variation of Current Collection with Time. 
All previous tests reported herein were for 
short durations, lasting at most a few minutes. 

A more extended test was conducted to determine 
the evolution that takes place with more pro- 
longed electron collection. Figure 9 shows that 
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Fig. 6. Comparison of a flat surface sample 
with samples having scribed patterns. 
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Fig. 7. The effect of sample temperature 
on electron collection. 
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Fig. 8. Data of Fig. 7 plotted to show the 
effect of conductor temperature at 
constant conductor potentials. 


the collected current does not decrease contin- 
uously with time, but shows occasional transients 
of large magnitude, often followed by gradual 
decreases. A visual inspection of the sample 
after the test showed that the polylmide insu- 
lator was darkened and beveled near the hole, as 
indicated in Fig. 10. The blackening and 
beveling of the sample supports the conclusion 
that, in time, the insulator is degraded by the 
energetic bombardment of electrons. This degra- 
dation should be included in any estimate of 
long-term collection effects in space. 



Fig. 9. Long-term variation of electron 
current collection with time at 
constant potential. 





(b) 


Fig. 10. Effect of current collection on 

insulator appearance, (a) untested 
(b) after 375 min. 

Variation of Insulator Thickness. Two 
insulator thicknesses were tested, 0.051 mm and 
0.127 mm. These were tested using four hole 
sizes: 0.35, 0.52, 1.0 and 2.0 mm diameter. 

For the samples tested there was no significant 
difference in current collection. It should be 
noted, though, that the insulator thickness was 
small relative to the hole diameter for all 
configurations tested. 
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Effect of Insulator Area. An insulator 
area effect had been found in an earlier inves- 
tigation.^ Tests were conducted to look for a 
similar effect in this investigation. To test 
for an area effect, three samples were con- 
structed: 2x2 cm, 5x5 cm, and 12 x 12 cm 

samples. These samples gave a surface area 
range of 4 cm^ to 144 cm^. The results of these 
tests are shown in Fig. 11. From this figure, 
it is seen that there is some area effect. The 
effect shown in Fig. 11, though, is smaller than 
that reported earlier by Kennerud.^ The nor- 
malized results shown lie within the extremes of 
his data which exhibit a much sharper rise in 
current collection with insulator area sur- 
rounding the hole. This difference is possibly 
due to the experimental differences mentioned 
earlier . 



Conductor Potentially 


Fig. 11. Comparison of samples with 

different insulator surface areas. 


Negative dias 


Normalization. Since ion collection 
involves different charge carriers than electron 
collection, a different normalization procedure 
should be expected. The current density of ions 
arriving at a planar boundary is given by the 
Bohm current density, H 



where m^ is the ion mass and all other parameters 
are as defined previously. The gas used for the 


hollow cathode (argon) is assumed for mj|^. As 
will be shown, a simple exchange of this current 
density for random electron current density did 
not result in adequate correlation of ion col- 
lection data. 


Effect of Pinhole Size. Four hole diameters 
were tested: 2.0, 2.95, 4.1 and 4.99 mm. The 

selection of such large hole sizes were based 
largely on the small ion currents collected. 

The data of these tests were normalized using 
both the Bohm current density, jBohm» 

Bohm current, Tgohm = ^hole^Bohm* see Figs. 

12(a) and 12(b). These results show that, 
unlike positive bias, hole area is an important 
parameter for ion current collection. Also, for 
the voltage range covered, the ion current 
collected varied nearly linearly with negative 
voltage. Preliminary tests at larger negative 
voltages than 1000 V have shown large, rapid, 
and variable current increases. The trends 
indicated in Figs. 12(a) and 12(b) should therefore 
not be extrapolated beyond the range shown. 




Conductor Potential, -V 


(a) 

Fig. 12. Ion current collection characteris- 
tics for a range of hole sizes, 
normalized by (a) ^Bohm* 
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(b) 

Fig. 12. Continued. 


Modeling 


For positive bias (electron collection), a 
crude model has been developed. This model is 
for the current collection process before any 
high-current mode change occurs as discussed 
earlier. It is assumed that any glow discharge 
effects are negligible in this lower voltage 
region, and that the collection process Involves 
only electrostatic collection from the plasma 
and a collection enhancement due to secondary 
electron emission on the insulator surface. The 
model is described by 


I = J. 


2 2 

1 + ^ + (1 _ ^)_V_ 

^ ^ 4 4^kT 


kole ^ ^ff 


(3) 


The term |l + \ + (1 “ \) -i^J is the 

collection term found by Parker and Whipple for 
collection by a planar probe surrounded by a 
conductor. The b in this term is a constant, 
found to be about 1,8 for the test environment 
used herein, V is the conductor potential, k is 


Boltzmann’s constant, Te Is the electron 
temperature, and jo is the random current density. 
The term Ahole is the area of exposed conductor. 
The term Agff (298/T)C(V) takes into account the 
collection involving secondary electron emission 
that gives the temperature dependence of the 
secondary electron emission, where T is the 
sample temperature (in ®K) , and 

^(V) = 2.1 + 4.3 X 10“^ V , (4) 


found experimentally. The term A^ff is the area 
of Insulation around the hole that contributes 
to the electron current collection, it also 
includes the effect of the secondary electron 
emission coefficient at 298®K. Agff is assumed 
as a function of only conductor potential. 

Using this model, the effective radius, r^ff, as 
measured from the hole edge, was calculated from 
the data and is shown in Fig. 13. The calcu- 
lated values of radius were then fitted to a 
simple parabolic variation with potential. The 
best least-square fit was obtained with 


This empirical fit was then used to calculate 
expected current/voltage curves. Typical cal- 
culated curves are shown in Fig. 14(a) and 14(b), 
together with the corresponding experimental 
data. 


o 



Fig. 13. Effective collection radius as a 
function of conductor potential. 

For negative bias (ion collection), the 
data in the voltage region investigated is in 
approximate agreement with the expression found 
by Parker and Whipple for a planar probe, 
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Conductor Potential, - V 


Fig. 15. Comparison of experimental ion 

current collection characteristics 
with planar probe theory (parameter 
b is an adjustable parameter of the 
probe theory) . , 

Concluding Remarks 


Fig. 14. Comparison of experimental and 
calculated electron current 
collection characteristics. 

2 2 

1 . _V K -V- /fi 

E “ 4 E (E +V) ^ 

o o o 

where Eq ** 8 kT^. This agreement is shown in 
Fig, 15 for different values of b. Note that 
the collection area for ions is simply the area 
of the hole. For the low voltage data presented 
for ion collection, this simple modeling tech- 
nique appears adequate. 


Bohm 


In this Investigation, two processes have 
been found to add to the current collection 
mechanism for positive bias. These were secon- 
dary electron emission by the Insulator surface 
surrounding the surface defect made apparent by 
the temperature dependence, and at higher vol- 
tages a glow discharge effect. The evidence 
supporting the secondary electron emission are 
the temperature effect and the effect of rough- 
ening the insulator surface. The glow discharge 
effect was easily seen at high currents, and the 
change of mode is probably related to some 
breakdown in the neutral background gas. 

Previous investigations where lower background 
pressures were used (<10“^ Torr) , did not report 
this effect. 


A model using secondary electron emission 
for the enhancement of current collection area 
has been developed. The empirical fit used for 
area enhancement gives a reasonable agreement to 
experimental results. Because of the empirical 
fit used, a change in insulator materials would 
require additional data for model verification. 

Negative bias experiments were done in a 
low voltage region where experimental results 
were consistent and agreed with planar probe 
theory. Although data were not presented at 
higher negative voltages, currents at these 
higher voltages were highly variable, with rapid 
transients that tended to damage instrumentation. 
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APPENDIX B. SECONDARY ELECTRON EMISSION FROM INSULATORS 

It has been believed that secondary electron emission from the 

insulating material surrounding the hole contributes to the electron 

1 2 

collection mechanism * (see Appendix A). With this in mind, modeling 
of the electron collection process will necessarily include some of the 
secondary emission properties. This appendix is a brief review of the 
properties that have been, or are felt to be, important to future 
modeling needs. 

Secondary Electron Emission Yields . Secondary emission yields 
consist of two components, true secondary electrons (emerging electrons 

3 

with energy less than 50 eV) and backseat tered primaries. The secondary 
emission yield is given by the ratio 

= number of emitted electrons /'n i 

number of incident electrons 

This total yield can also be written as 

a = 6 + n , (B-2) 


where 6 is the yield of true secondaries 

5 _ number of true secondary electrons emitted 
number of incident electrons 


(B-3) 


and n is the backscattered yield 

= number of backscattered electrons . , . 

^ number of incident electrons * ^ ^ 

The total yield is dependent on the incident electron energy, but 6 and 
f) exhibit different dependencies. 
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In the present investigation, the insulating materials of interest 

are polyimide and polytetraf luoroethylene (teflon). Studies of emission 

4 

yields (a) of these materials have been done by Willis and Skinner. 

The secondary emission yield (6) of polyimide has been investigated by 
Leung, Tueling and Schnauss.^ The emission yield (a) for polyimide from 
Willis and Skinner (angle of incidence of primary electrons not given) 
and the secondary yield (o) from Leung, Tueling and Schnauss (for normal 
incidence) are shown in Fig. B-1. 

The emission yields from Willis and Skinner have been fit to the 

f orm^ 


a(E) = K 



-bE 

e 


) 


(B-5) 


where E is the energy of the primary electron in eV and K, a and b are 

constants dependent on the material. For polyimide, K = 3.5, 

-3 -2 

a = 2 X 10 and b = 1.5 x 10 , and for polytetrafluorethylene (teflon) 

-3 -3 

K = 5.8, a = 1 X 10 and b = 5 x 10 . The secondary yield from Leung, 

Tueling and Schnauss is given by 

6(E) = exp ~ - \/2 ~ (B-6) 

Eli E 

max 

where Is the primary energy where 6 is maximum (250 eV for polyimide) 

and E^^ is the primary energy where 6=1 (the second crossover point, 

E^^ = 1000 eV for polyimide, see Fig. B-1). All energies are in eV. 

Please note that since the primary electron reflection coefficient n 
can vary between 0 and 1, a should at all times be greater than 6. 

The amount by which this is not true for the data described (Fig. B-1) 
indicates the probable error in this data. 



Yield 
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Fig. B-1. (a) Total yield from Willis and Skinner^ and (b) secondary 

electron yield from Leung, fueling and Schnauss^ as a 
function of incident electron energy. 
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Angular Dependence of Secondary Emission Yields , Secondary emission 

yields are dependent on the incidence angle of the primary electrons. 

9 

This dependence is given by 

5 = 6 (B-9) 

0 o 

where 0 is the angle of incidence of electrons with respect to the 
surface normal, 6^ is the secondary emission yield at normal incidence, 
is the secondary emission yield for electrons incident at an angle 0, 
and p is a constant of the material independent of 0. The constant p 
is related to the mean depth of penetration of primary electrons and 
the absorption coefficient of secondary electrons. It has been deter- 

g 

mined empirically to be approximately 2 for most polymers. 


Energy Distribution of Secondary Electrons . A typical energy 
distribution is shown in Fig. B-2. Although this distribution is for 
a metal, many insulators have been found to have similar distributions. 
The energy distributions have been found to be independent of primary 
electron energy range of 20 to 1000 eV.^^ 


10 


Angular Distribution of Secondary Electrons . The secondary electrons 

12 

are emitted from the surface with approximately a cosine distribution 
which is nearly independent of the angle of incidence of the primary 
electrons. 


Effect of Temperature . The secondary yield has been found to have 

13 11 

a temperature dependence for insulators, although not for metals. 

The reason for this difference appears to be in the energy loss mechanism 
of the secondary electrons in the material. In metals, the secondaries 
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can collide with free and bound electrons involving large and small 
amounts of energy transfer . In insulators, because of the large gap 
between the valence and conduction bands, any energy loss mechanism of 
the secondaries must depend on electron collisions with lattice defects 
and/or electron-phonon interactions for energies less than the band 
gap.^^ For ionic crystals, a temperature dependence of 



2 \ 1/2 

T + 1 

e^^/kT2 - 1 

-^1 

e^^/kT^ - 1 


(B-8) 


has been derived, and is in good agreement with experimental results. 
In the expression, h is Planck’s constant, k is Boltzmann’s constant, 

T is the insulator temperature and v is the frequency of optical 
vibrations. For kT >> hv, this expression (Eq. B-8) becomes 



which has been mentioned in a previous report. 

Non-polar insulators have a temperature dependence of 


(B-9) 


max/ 


f 


1/2 


e^^/kT -1 
o \ 2 


+ 1 V dv 


‘'max , 


(B-10) 


/ 


le^^/kT-1 


+ v^dv 


For kT >> hv, this expression reduces to Eq. B-9. 
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Effect of Surface Condition . It has been found that when a surface 
is roughened, the secondary yield decreases. The explanation for this 
is that when a surface is roughened, microscopic Faraday spaces (spaces 
without fields and from which no secondary electrons can be emitted) 

are formed. Thus decreasing the total area from which emission may 

16 2 
occur. This effect has also been mentioned in a previous report. 
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